Strongly Correlated Superconductivity and Pseudogap Phase near a multi-band Mott 
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Near a Mott transition, strong electron correlations may enhance Cooper pairing. This is demon- 
strated in the Dynamical Mean Field Theory solution of a twofold-orbital degenerate Hubbard model 
with an inverted on-site Hund rule exchange, favoring local spin-singlet configurations. Close to the 
Mott insulator, which here is a local version of a valence bond insulator, a pseudogap non-Fermi- 
liquid metal, a superconductor, and a normal metal appear, in striking similarity with the physics 
of cuprates. The strongly correlated s-wave superconducting state has a larger Drude weight than 
the corresponding normal state. The role of the impurity Kondo problem is underscored. 

PACS numbers: 74.20.Mn, 71.27.+a, 71.30.+h, 71.10.Hf 
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How superconductivity emerges so spectacularly out of 
a weakly doped Mott insulator is one of the fascinating 
and still controversial aspects of high Tc cuprate super- 
conductors. A novel strongly correlated superconductiv- 
ity (SCS) scenario has been recently proposed Q which 
deals with an ultimately related basic issue, namely the 
conditions under which Cooper pairing can be enhanced, 
rather than depressed, by strong electron repulsion. The 
key of the SCS proposal is the presence of pairing in- 
teraction term J, which is weak but is not suppressed 
by the strong short-range repulsion U. This is real- 
ized for an attraction which involves mainly spin and 
orbital degrees of freedom, which are not frozen near a 
Mott metal-insulator transition (MIT) . In addition, close 
to the MIT, correlations slow down electron motion so 
much that the effective quasiparticle bandwidth becomes 
extremely small Wqp W, W being its uncorrelated 
value. In these conditions the pairing attraction can 
eventually equal the quasiparticle bandwidth J ^ Wqp. 
That drives the system to an intermediate-strong cou- 
pling superconducting regime where the maximum su- 
perconducting gap A ~ J for given value of J is reached, 
as opposed to the much smaller uncorrelated BCS value 
^BCS ~ W exp{—W/J). A first theoretical realization 
of SCS was demonstrated in Ref. 0, 13| by a Dynamical 
Mean Field Theory (DMFT)[i| solution of a three-fold 
degenerate model for tetravalent fuUerides. It was shown 
that a narrow SCS region arises next the MIT, and is 
indeed characterized by a superconducting gap three or- 
ders of magnitude larger that the value one would obtain 
for the same attraction in the absence of correlation. 

DMFT maps a lattice model onto an Anderson impu- 



rity (AT) model subject to a self-consistency condition. 
Within this mapping SCS emerges through a compe- 
tition between Kondo screening of the impurity (leading 
to a normal Fermi liquid) and an intra-impurity mecha- 
nism which forms a local non-degenerate singlet (a kind 
of local resonating valence bond state 01 )• This requires 
orbital degeneracy and inversion of Hund's rules, both 
ingredients present in the model of Ref. [IJ. More re- 
cently a simpler AI with only twofold orbital degeneracy 
was shown by Wilson Numerical Renormalization Group 
(NRG) to display anomalous properties (5|], suggesting its 
lattice generalization as a new candidate for a SCS. 

In this Letter we present a detailed DMFT analysis 
that confirms this expectation, exploiting the lower de- 
generacy for a wider and more revealing study. We con- 
sider an infinite coordination Bethe lattice and solve the 
AI by exact diagonalization. The model reads 
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where cj creates an electron at site i in orbital a = 

1, 2 with spin a, while Ui — J^a a "^I aa'^i aa ^be on-site 
occupation number. The on-site exchange is 
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where T^^^ = 1/2 Ea.fc 4,a^^ (^")afc c,_6<x are the 
pseudo-spin 1/2 operators, and f" (a = x,y,z) the 
Pauli matrices. The electronic states of the isolated site 
with n electrons are labeled by total spin and pseudo- 
spin, S and T and their z-components, and have energies 
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E{n,S,S„T,T,) = Un'^/2 - 2J[T{T + 1) - T^]. For 
n — 2 the configurations allowed by Pauli principle are 
a spin triplet orbital singlet (5 = 1 and T = 0) and a 
spin singlet orbital triplet [S — and T — 1), split by 
Hj. If J < 0, standard Hund's rules, the spin triplet 
has the lowest energy. Here we consider the less com- 
mon case of J > 0, where the lowest energy configura- 
tion has S ~ 0, T = 1 and = 0. This inversion of 
Hund's rules may for instance mimic a dynamical e® E 
Jahn-Teller effect HQ- Here it just represents a generic 
mechanism for on-site spin-singlet pairing. Indeed, for 
[/ = and J ^ W, the ground state of Q is an s- 
wave BCS superconductor with pairs condensed in the 
5 = 0, T = 1 and — Q channel. The energy gap is 
A ~ VFexp(— 1/A), where A — 2JNp is the dimensionless 
superconducting coupling, and Np the density of states 
(DOS) at the Fermi energy per spin and orbital. In this 
regime, a finite U <^ W introduces a "Coulomb" pseu- 
dopotential fi^ ~ U Np that opposes superconductivity, 
eventually suppressed for U > 2 J in favor of a normal 
metal ground state. For larger U > Uc ^ W, the model 
undergoes a MIT for all integer fillings (n) = 1,2,3. 
We could generally expect the superconducting gap to 
decrease monotonically as a function of U, the supercon- 
ductor either turning directly into a Mott insulator (for 
A ~ 1), or (when A ^ 1) first into a metal, and then into 
the insulator. This is indeed what we find for (n) = 1 
(equivalent to (n) = 3 by particle-hole symmetry) and A 
ranging from 0.1 to 0.6. The picture is however richer in 
the half-filled (n) = 2 case, as reported in Fig. ^ where 
the standard superconductor evolves continuously to the 
SCS regime If A is large (e.g.,J/M^ = 0.15), A decreases 
with increasing U, as in a BCS superconductor, until a 
weakly first order transition to our local-RVB insulator 
occurs 6]. For smaller A however, the dependence of A 
on U is non-monotonic, and the initial drop for small U 
is followed by a rise. Finally, for the smallest values of A, 
the weak-C/ superconductor first turns into a metal but, 
just before the MIT, reverts back to a superconductor, 
as in SCS of Ref. 0. Here too the superconducting gap 
reaches values much larger than those attained for the 
same pairing attraction A but for U = (compare the 
inset and the right-hand side of Fig.^. 

The appearance of SCS in |^ for (n) = 2 supports the 
AI analysis 0. Indeed, as U increases, the quasiparticle 
bandwidth, which is the AI effective Kondo temperature 
Tk , is gradually suppressed by the dropping quasiparticle 
residue Z <C 1, Tk ^ ZW. At the same time, we expect 
J to remain unrenormalized . Hj splits multiplets 
with same n, and that can happen without opposing f/, 
which just freezes charge fluctuations regardless of either 
spin or orbital degrees of freedom. As Z decreases, the 
AI enters the critical region around its unstable fixed- 
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point 13 at a critical tI^^ - J. For Tk > T'^' , the AI 
is in the Kondo screened regime, while for Tk < T^' 



FIG. 1: Superconducting gap for (n) = 2 as a function of 
Hubbard repulsion for fixed coupling J= X/{2Nf) ■ Increas- 
ing repulsion spoils superconductivity at large coupling. At 
weak coupling superconductivity is instead strongly enhanced 
close to the Mott transition. The inset shows the weak cou- 
pling regime on an expanded scale, showing a much smaller 
gap at small U compared with SCS at large U. 



an intra-impurity singlet forms thanks to the inverted 
exchange J. At the fixed point the two effects balance 
exactly, leading to a residual entropy In \/2. The remain- 
ing impurity degrees of freedom are quenched at a larger 

temperature T+ ~ max{TK, J). r_ - \Tk ~ T^k^/^k'' 
measures the deviation from the fixed point behavior. 
It was argued in Ref. that DMFT self-consistency 
should turn this AI instability into a true bulk one, most 
likely implying superconductivity. This is now confirmed 
by our DMFT analysis. S-wave superconductivity opens 
up a new screening route which freezes the residual en- 
tropy, otherwise quenched only below r_'5'|. This sug- 
gests that (a) the energy scale which controls supercon- 
ductivity should be related to r+ — r_; (b) the onset 
of superconductivity should be accompanied by a kinetic 
energy gain at low frequency, as Kondo screening im- 
plies. A kinetic energy (or, more accurately, band en- 
ergy) gain may be regarded as a signal of SCS (as op- 
posed to BCS where kinetic energy rises), and that is 
reflected by the behavior of the Drude weight. In Fig. [3 
we show for J — 0.05 W the Drude weight of the stable 
superconducting solution for 0.7 < U /W < Uc (here the 
Drude weight is the strength of the superfluid peak), in 
comparison with that of the unstable metallic solution 
(obtained disallowing superconductivity), taken to rep- 
resent the normal phase. Fig.|21shows a large increase of 
Drude weight with superconductivity, an occurrence also 
predicted 0, 01 and actually observed in cupratesfl^. 
Here the unstable metal Drude weight actually appears 
to vanish at J7* < Uc- As we shall show, this reflects the 
opening of a pseudogap at the chemical potential before 
the MIT. 

Additional physical insight is gained by analyzing the 
DOS in the strongly correlated region U/W > 0.7. 
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FIG. 2: Drude weight as function of [/ at J = 0.05W. 
SCS points: lowest energy superconducting solution; unstable 
metal points : metastable solution with superconducting or- 
der parameter forced equal to zero. Note the very large Drude 
weight increase in the SCS pocket (magnified in the inset). 
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FIG. 3: Single-particle DOS at J = 0.051V of the metallic 
solution for diff'erent U's close to the MIT (superconductivity 
between U = 0.7 and Uc — 0.88 is not allowed). Note the 
pseudogap within the low-energy quasiparticle peak coexist- 
ing with the high-energy Hubbard bands. 



Within DMFT, the DOS in this region of parameters dis- 
plays three features: Two high-energy Hubbard bands 
and a low-energy feature associated to quasiparticles. 
Fig. O displays the evolution of the DOS of the normal 
state solution in this regime, and shows that a pseudogap 
opens within the low-energy peak. A similar pseudogap 
appeared also in the NRG solution of the AI of Ref . [3 , 
where it was argued that the low-energy DOS around the 
unstable fixed point is well described by 

where po is the non-interacting DOS, and plus/minus 
refers to the Kondo screened/unscreened phase. We find 
that (PI fits well also the low-frequency behavior of the 
DMFT results 's'l. Eq. © impHes that at the chemical 
potential p(0) = po in the Kondo screened regime, com- 
patible with a Fermi-liquid with XmS](w) ~ cj^. On the 
contrary, the DOS values at r_ = 0, p(0) = po/2, and in 
the pseudogap phase, p{0) = 0, imply a singular behavior 
of the self-energy, i.e., a non- Fermi liquid. The best-fit 
values of T+ and for J = 0.05 W are shown in Fig.^ 
and compared with the superconducting gap A. The 
maximum of A almost coincides with the vanishing of r_ , 
which corresponds to the AI unstable fixed point. Fig. 0] 
suggests a scenario which shares similarities to the quan- 
tum critical point 0| and the gauge theory-slave boson fioj 
pictures of cuprates. For temperatures T above Tc (pre- 
sumably of order A), but below T_, the normal phase is 
a Fermi-liquid for U < U^:, and a non Fermi-liquid pseu- 
dogap phase for Uc > U > U^. For T. < T < T+, both 
the narrow Fermi-liquid quasiparticle peak and the pseu- 
dogap should be washed out, leaving a broader resonance 
reflecting the properties of the non Fermi-liquid AI un- 
stable fixed point. The resonance will eventually disap- 



pear above T^. Close to the MIT, the unstable metallic 
solution and the stable superconductor have almost the 
same energy and very similar DOS, but for the presence 
of a very small superconducting gap. From this point of 
view, the pseudogapped unstable phase plays a role simi- 
lar to the staggered flux phase within the SU(2) invariant 
slave-boson description of the t-J model, and may there- 
fore be thought as a phase with broken symmetry in the 
particle-hole or in the particle-particle channels, where 
the full symmetry is restored by fluctuations^] . We note 
here that (a) superconductivity is not an accidental route 
which the lattice system takes to rid itself of competition 
among other phases, but is one of the pre-determined 
instabilities of the AI; (b) attempts to uncover the fixed 
point by suppressing superconductivity would likely spoil 
the fixed point or unveil other instabilities of the lattice 
model. 

The analogy with high Tc cuprates becomes even more 
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FIG. 4: Behavior of the relevant energy scales (defined in the 
text) close to the MIT for J/W = 0.05 as function of U/W. 
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FIG. 5; Phase diagram as function of U/W and doping 
S — n — 2sLtJ = 0.05 W . The think vertical line marks the 
singlet Mott insulator. The inset shows, for U = 0.92Vy, the 
superconducting gap A divided by a factor 10"'^ and Drude 
weight D (normalized to the non-interacting value) as func- 
tions of doping 5. 



suggestive when we analyze the phase diagram away from 
half-fiUing, and foUow the fate of superconductivity upon 
doping. SCS extends (Fig. [S)) into a superconducting 
pocket away from {n) — 2. When we dope the singlet 
Mott insulator just above the MIT, a pseudogap metallic 
phase is encountered first, followed by a narrow supercon- 
ducting region which finally gives way to a normal metal 
(inset of Fig.|3J), in remarkable similarity to the phase di- 
agram of high Tc cuprates. Striking is also the similarity 
in the Drude weight (inset), which is zero in the Mott in- 
sulator at half-filling and increases almost linearly upon 
doping. The phase diagram of Fig. |31 again agrees with 
the AI analysis 0, which shows that the unstable fixed 
point still exists away from half-filling. 

In the alkali fulleride models, that provide our prime 
examples of SCS, the physics which leads to J > and 
spin pairin g is Jahn Teller, i.e. a conventional phonon 
effect, [ilia, im Future experiments to test the SCS sce- 
nario by comparing Drude weights in presence and in 



absence of superconductivity should be actively consid- 
ered in trivalent fuUerides such as K'^NH^Cfin. jl2j These 
materials are of extraordinary importance, for they show 
a similar "universal" contiguity between superconductiv- 
ity and Mott insulating behavior as the high Tc cuprates 
do, albeit with different physical ingredients. When deal- 
ing with spectroscopies, one should keep in mind that 
there are now two low-energy scales, T+ and That 
makes, e. g.. p hotoemission data such as the recent ones 
on K3C60I13 especially difficult to analyze, even ignor- 
ing the heavy vibronic structures. In the SCS region, 
coherent band-like dispersion should occur with a re- 
duced bandwidth of order 2T_|_, a scale which remains 
finite across the MIT, even though the system is strongly 



correlated. The true quasiparticle peak will instead be- 
come extremely narrow, of order T_, or even disappear 
in the pseudogap regime. We believe that the datafisj 
are actually compatible with the SCS scenario. The un- 
derlying Kondo physics offers a new intriguing insight in 
this superconductivity problem. It seems in fact that the 
onset of bulk coherence close to a MIT may involve a 
hierarchy of energy scales, T- < T+ and A. The latter 
is the most sought-after, as it marks the onset of lattice 
long-range order. On the contrary, T+ and r_ have their 
clearest meaning in the AI, which displays a two-stage 
quenching of the impurity entropy. 
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